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Abstract 



These final results on e~^e~ — > W^W~ production cross-section measurements 
at LEP2 use data collected by the DELPHI detector at centre-of-mass energies 
up to 209 GeV. Measurements of total cross-sections, W angular differential 
distributions and decay branching fractions, and the value of the CKM element 
\Vcs\ are compared to the expectations of the Standard Model. 
These results supersede all values previously published by DELPHI. 
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1 Introduction 

The cross-section for the doubly resonant production of W bosons has been measured 
with the data sample collected with the DELPHI detector during the high-energy opera- 
tion phase of the LEP e^e~ collider (LEP2), at centre-of-mass energies from 161 GeV up 
to 209 GeV. Measurements of total cross-sections, W angular differential distributions and 
decay branching fractions, and the value of the CKM element \Vcs\ are presented. They 
are compared to the Standard Model including most recent theoretical predictions [1]. 

The cross-sections determined in these analyses correspond to W pair production, 
defined through the three doubly resonant tree-level diagrams ("CC03 diagrams" [2] in 
the following) involving s-channel 7 and Z exchange and t-channel v exchange, as shown 
in Figure 1. Depending on the decay mode of each W, fully hadronic, mixed hadronic- 
leptonic ( "semi-leptonic" ) or fully leptonic final states arc obtained. 

This paper is organised as follows: after a brief description of the DELPHI detector 
in Section 2, a summary of the analysed data and simulation samples is provided in 
Section 3. Track selection and particle identification are briefiy illustrated in Section 4 and 
in Section 5 the selection of WW events is described and the performance of the analysis 
reported; Section 6 is devoted to the discussion of the systematic error assessment. In 
Sections 7 and 8, results in terms of differential cross-sections, total cross-section and W 
branching ratios are presented. In Section 9 the results are compared to the theoretical 
predictions and conclusions follow in Section 10. 

2 The DELPHI detector in the LEP2 phase 

The DELPHI detector configuration for LEPl running was described elsewhere [3]. 
For operation at LEP2, changes were made to the subdctcctors, the trigger system [4], 
the run control and the algorithms used in the ofllinc reconstruction of charged particles, 
which improved the performance compared to LEPl. 

The major changes were the extension of the Vertex Detector (VD) and the inclusion 
of the Very Forward Tracker (VFT) [5], which enlarged the coverage of the silicon tracker 
out to 11° < 6 < 169°^. Also the Inner Detector, both the Jet Chamber and Trigger 
Layers, were extended to cover the polar angle region 15° < 9 < 165°. Together with 
improved tracking algorithms, alignment and calibration procedures optimised for LEP2, 
these changes led to an improved track reconstruction efficiency in the forward regions. 

Changes were made to the electronics of the trigger and timing system which improved 
the stability of the running during data taking. The trigger conditions were optimised for 
LEP2 running, to give high efficiency for Standard Model two-fermion and four-fermion 
processes and also give sensitivity to events which may be signatures of new physics. 
In addition, improvements were made to the operation of the detector during the LEP 
cycle, to prepare the detector for data taking at the very start of stable collisions of the 
beams, and to respond to adverse background conditions from LEP were they to arise. 
These changes led to an overall improvement of ~ 10% in the efficiency for collecting the 
dehvercd luminosity, from ~ 85% at the start of LEP2 in 1995 to ~ 95% at the end in 
year 2000. 

^The DELPHI coordinate system has the 2-axis coUinear with the incoming electron beam. indicates the polar angle 
with respect to the 2-axis, Rep indicates the plane perpendicular to the z axis. 
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During the operation of the detector in year 2000 one of the sectors representing 
1/12 of the acceptance of the central tracking device, the Time Projection Chamber 
(TPC), failed. This problem affects about a quarter of the data collected in that year. 
Nevertheless, the redundancy of the tracking system meant that charged particles passing 
through the sector could still be reconstructed from signals in other tracking detectors. 
A modified track reconstruction algorithm was used in this sector, which included space 
points reconstructed in the Barrel RICH detector, helpful in the determination of the 
polar angle of charged particles. As a result, the track reconstruction efficiency was only 
slightly reduced in the region covered by the broken sector. The impact of the failure of 
this part of the detector on the analyses is discussed further in Section 6.1. 

3 Data and simulation samples 

A summary of the data samples used for the WW cross-section measurement is re- 
ported in Table 1, where the luminosity- weighted centre-of-mass energies and the amount 
of data analysed at each energy are shown. The luminosity is determined from Bhabha 
scattering measurements making use of the very forward electromagnetic calorimetry [6] . 
The total integrated luminosity for the LEP2 period corresponds approximately to 
670 pb~^. The luminosities used for the different selections correspond to those data 
for which all elements of the detector essential to each specific analysis were fully func- 
tional; tighter requirements on the detectors used for Icpton identification were applied 
for the data samples used in the semi-leptonic and fully leptonic channel analyses. The 
luminosity in year 2000 was delivered in a continuum of energies, thus data taken during 
this period arc divided into two centre-of-mass energy ranges, above and below 205.5 GeV, 
referred to as 205 GeV and 207 GeV in the following. 

All the data taken from the year 1997 onwards have been reprocessed with an updated 
version of the DELPHI reconstruction software, using a consistent treatment for all the 
samples. Larger simulation samples were realised with more up-to-date Monte Carlo pro- 
grams, interfaced to the full DELPHI simulation program DELSIM [3] and reconstructed 
with the same reconstruction program as the real data. 

The cross-section analyses on these data arc updated with respect to the previously 
pubhshed ones [7] and supersede them. The data taken in year 1996 have not been re- 
analysed, because possible improvements were negligible compared to the large statistical 
errors of those measurements; these results correspond to the publications in [9], with a 
revised determination of the luminosity [8]. 

Four-fermion simulation samples were produced with the WPHACT [10] generator, 
interfaced with the PYTHIA 6.156 [11] hadronisation model. In order to perform checks 
on hadronisation effects, the same generator was also interfaced to the ARIADNE [12] 
and HERWIG [13] hadronisation models. 

The generation was complemented with two-photon collision generators BDK [14], 
BDKRC [15] and PYTHIA. The most recent 0{a) electroweak corrections, via the so- 
called Leading Pole Approximation (LPA), were included in our generation of the signal 
via weights given by the YFSWW program [16], according to the scheme described in [17]. 

The selection efficiencies were defined with respect to the CC03 diagrams only, by 
reweighting each event to the CC03 contribution according to the ratio of the squared 
matrix elements computed with these diagrams only and with the full set of diagrams. 

The simulation of two-fermion background processes was realised with the KK2f [18] 
and KoralZ [19] generators interfaced to PYTHIA, ARIADNEand HERWIG, for 
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Year 


£- weighted ^/s (GeV) 


Hadronic JC (pb ^) 


Leptonic £ (pb ^) 


1996 


161.31 


10.07 


10.07 




172.14 


10.12 


10.12 


1997 


182.65 


52.51 


51.63 


1998 


188.63 


154.35 


153.81 


1999 


191.58 


25.16 


24.51 




195.51 


76.08 


71.99 




199.51 


82.79 


81.82 




201.64 


40.31 


39.70 


2000 


204.81 


82.63 


74.93 




206.55 


135.82 


123.66 



Table 1: Integrated luminosity- weighted centre-of-mass energies and luminosities in the 
LEP2 data taking period. The hadronic luminosity is used for the four-quark channel, 
the leptonic one for the semi-leptonic and fully leptonic channels. 

which the fragmentation parameters have been tuned at the Z-resonance [20], and the 
BHWIDE [21] generator for Bhabha events. 

4 Charged particle selection and lepton identifica- 
tion 

To improve the event reconstruction and reject contributions from either cosmic ray 
events, beam-gas interactions or noise from electronics, the reconstructed charged parti- 
cles were required to fulfil the following criteria: 

- momentum 0.1 GeV/c <p< 1.5 • Pbeam', 

- momentum error Ap/p< 1; 

- i?0 impact parameter < 4 cm; 

- z impact parameter < 4 cm/sin^. 

Tracks seen by only the central tracking devices (Vertex Detector (VD) and Inner Detector 
(ID)) were rejected if no z coordinate measurement was available. 

Neutral clusters of energy in the barrel (HPC), forward (FEMC) and very forward 
(STIC) electromagnetic calorimeters were required to have an energy of at least 300, 
400 and 300 MeV respectively. In addition, off-momentum electrons in the STIC were 
rejected by a cut on the polar angle of a shower at 3°. Noise from the hadron calorimeter 
(HCAL) was reduced by rejecting showers which had activity in only one layer and failed 
a selection based on energy, number of hits and longitudinal position. 

Muon identification was performed by extrapolating tracks through the entire detector, 
and associating nearby HCAL energy deposits and muon chamber hits to the tracks. The 
errors on track parameters, on the energy deposits and on the position of chamber hits 
were taken into account when making this association. Particles were identified as muons 
if there was at least one muon chamber hit associated to a track or if the size and 
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longitudinal profile of the HCAL energy deposits associated to a track were consistent 
with a minimum ionizing particle. 

Electron identification was based on the reconstructed showers in the electromagnetic 
calorimeters associated to charged particle tracks. The barrel photons in a conic region 
close to the direction of the track were reassociated to the candidate electron, also ac- 
counting for the bending caused by the magnetic field. A particle was identified as an 
electron if the associated calorimetric energy was larger than 50% of the momentum re- 
constructed in magnetic field. For tracks below 30 GeV/c, the identification algorithm 
was based on a combination of the energy-momentum ratio, the measurement of the en- 
ergy loss in the TPC, the matching of the track to the shower in the z and directions 
and on the shower profile. FEMC electromagnetic deposits close in space were clustered 
together and the association with reconstructed tracks was used for electron/photon dis- 
crimination. Care was taken in excluding those tracks which were likely to come from 
the development of a shower outside the calorimeter. In addition, the ratio between the 
electromagnetic energy and the total energy (electromagnetic and hadronic) was required 
to be above 90%. For energetic and isolated tracks in the regions not covered by elec- 
tromagnetic calorimetry, the measurement of the energy loss in the TPC was used. The 
estimate of electron energy was obtained from a combination of the track momentum and 
the calorimetric deposition. 

5 Event selection and partial cross-sections 

5.1 Fully hadronic final state 

5.1.1 Selection of fully hadronic final state events 

For the selection of W+W^ qqqq, the charged and neutral particles in each event 
were forced into a four-jet configuration with the DURHAM algorithm [22]. A preselection 
of events was performed to select well reconstructed hadronic events without missing 
energy and momentum. The following conditions were imposed: 

• reconstructed effective centre-of-mass energy [23] a/s' > 65% of the nominal centre- 
of-mass energy; 

• total and transverse energy for charged particles, Ech = and Et = J2i \ pi \, 
where p] is the momentum component of the particle i perpendicular to the beam 
axis, each > 20% of nominal centre-of-mass energy; 

• total particle multiplicity > 3 for each jet; 

• Vcut > 0.0006 for the migration from 4 to 3 jets when clustering with the DURHAM 

algorithm; 

• convergence of a four- constraint (4C) fit of the measured jet energies and directions 
imposing four-momentum conservation. 

A feed-forward neural network was then used to improve the rejection of two-fermion 
(mainly Z/7 — qq(g)) and four-fermion background (mainly ZZ — > qqqq, qqT+r~). The 
network, based on the JETNET package [24] , uses the standard back- propagation algo- 
rithm and consists of three layers with 13 input nodes, 7 hidden nodes and one output 
node. 

The following jet and event observables were chosen as input variables, taking into 
account previous neural network studies [25] to optimise input variables for the WW and 
two-fermion separation: 
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• the difference between tlie maximum and minimum jet energies after tlie 4C fit; 

• tfie minimum angle between two jets after tfie 4C fit; 

• tlie value of ycut from the DURHAM algorithm for the migration of 4 jets into 3 jets; 

• the minimum particle multiplicity of any jet; 

• the reconstructed effective ccntre-of-mass energy 

• the maximum probability, amongst the three possible jet pairings, for a 6C fit (im- 
posing 4-momentum conservation and the invariant mass of each jet pair equal to 
the W-mass); 

• the thrust; 

• the sphericity; 

• the transverse energy; 

• the sum of the cubes of the magnitudes of the momenta of the 7 highest momentum 
particles Yli=i 

• the minimum jet broadening i?mm [22]; 

• the Fox- Wolfram-moment H3 [26]; 

• the Fox- Wolf ram-moment H4. 

The neural network was trained separately at each energy. Each training was per- 
formed with 2500 signal events and 2500 Z/7 — > qq background events. Afterwards the 
network output was calculated for other independent samples of simulated four-fcrmion, 
Z/7 and 77 events, and for the data. Figure 2 shows distributions of the neural network 
output for data and simulated events at 189 and 207 GeV. 

Events were selected by applying a cut on the NN output variable, chosen to minimise 
the total error on the measured cross-section including the systematic uncertainty on the 
two-fermion background with its correlation at all centre-of-mass energies (see Section 6). 

5.1.2 Results for fully hadronic final state 

The efficiency and background contamination for the hadronic event selections were 
evaluated independently at the different centre-of-mass energies. The selection perfor- 
mance at y/s — 200 GeV and the total number of events selected in each data sample are 
reported in Table 2. The efficiencies varied by no more than 4% over the different energy 
points above 172 GeV. The background is dominated by qq(7), representing 70-75% of 
the contamination, decreasing with energy. 

At each energy point the cross-section for fully hadronic events was obtained from a 
binned maximum likelihood fit to the distribution of the NN output variable above the 
cut value, assuming Poissonian probability density functions for the number of events. 
The probability is calculated on the basis of the efficiency for being reconstructed in a 
given bin of the NN output and the expected background in each bin and is a function of 
the partial cross-section to be measured. For this fit the contamination from other WW 
channels, with its value fixed to the SM prediction, was considered as a background. 

The results for = cr^mX BR (WW qqqq), where BR(WW qqqq) is the 

probability for the W-pair to give a purely hadronic final state, are reported in Table 3. 
Systematic uncertainties were determined as detailed in Section 6. 
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efficiencies for selected channels 


channel 


JJJJ 








qqqq 


0.797 


< 10-4 


< 10-4 


0.012 


qqev 


0.004 


0.677 


0.004 


0.114 


qqjiv 


0.002 


0.001 


0.852 


0.043 


qqru 


0.016 


0.032 


0.026 


0.581 


background (pb) 


1.21 


0.232 


0.075 


0.371 


(GeV) 


Selected events 


161 


15 




12 




172 


65 


14 


17 


14 


183 


345 


94 


118 


123 


189 


1042 


269 


336 


339 


192 


187 


42 


53 


58 


196 


532 


151 


166 


164 


200 


614 


162 


190 


208 


202 


317 


94 


89 


83 


205 


657 


169 


153 


174 


207 


999 


214 


259 


289 


All 


4773 


4054 



Table 2: Data for the cross-section measurement of hadronic and semi-leptonic final 
states. The efficiency matrix and the backgrounds are the ones at 200 GeV. The back- 
grounds include two-fermion and non-CC03 four-fermion contributions. The upper limits 
on the efficiencies are at the 95% C.L.. 

5.2 Semi-leptonic final state 

5.2.1 Selection of semi-leptonic final state events 

Events in which one of the W bosons decays into a lepton and a neutrino and 
the other one into quarks are characterised by two or more hadronic jets, one iso- 
lated lepton (coming either directly from the W decay or from the cascade decay 
W — > eVfiVrfT- oi /iz/^z/^z/,-) or a low multiplicity jet due to a hadronic r decay, 

and missing momentum resulting from the neutrino (s). The major background comes 
from gg(7) production and from four-fermion final states containing two quarks and two 
leptons of the same flavour. 

Events were first required to pass a general hadronic preselection: 

• at least 5 charged particles; 

• energy of charged particles at least 10% of total centre-of-mass energy; 

• yjEMFf + EMF^ < 0.9 x Eheam, where EMFf^b identify the total energy deposited 
in electromagnetic calorimeters in the forward and backward directions, defined as 
two 20° cones around the beam axes. 
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(GcV) 


^ww^ = c^ww X BR(WW qqqq) (pb) 


161 


1.53toil (stat) ±0.13 (syst) 


172 


4.65i[J:i (stat) ± 0.18 (syst) 


183 


7.23 ± 0.45 (stat) ± 0.09 (syst) 


189 


7.38 ± 0.27 (stat) ± 0.09 (syst) 


192 


7.78 ± 0.68 (stat) ± 0.10 (syst) 


196 


7.69 ± 0.39 (stat) ± 0.10 (syst) 


200 


7.73 ± 0.37 (stat) ± 0.10 (syst) 


202 


7.83 ± 0.54 (stat) ± 0.10 (syst) 


205 


8.26 ± 0.38 (stat) ± 0.10 (syst) 


207 


7.59 ± 0.29 (stat) ± 0.10 (syst) 



Table 3: Measured hadronic cross-sections. 



A search for leptons was then made. Of the identified electrons with an energy greater 
than 5 GeV, the one with the highest value of (energy x ^j^o) ^ was considered to be the 
electron candidate. This candidate was required to have an energy of at least 15 GeV. Of 
the identified muons with a momentum greater than 5 GeV/c, the one with the highest 
value of (momentum x 9iso) was considered to be the muon candidate. This candidate 
was required to have a momentum of at least 15 GeV/c. The event was then clustered 
into jets using the LUCLUS [11] algorithm with a djoin value of 6.5 GeV/c. The resulting 
jets were trimmed by removing particles at an angle greater than 20° to the highest 
energy particle in the jet. Of these trimmed jets, the one with the smallest momentum- 
weighted spread ^ was taken to be the tau candidate. Particles with momenta smaller 
than 1 GeV/c were removed from the candidate jet if they were at an angle greater than 
8° to the jet axis. At the end of the procedure, this jet was required to still contain at 
least one charged particle. 

An event could have up to three lepton candidates, one of each flavour. 

For each lepton candidate, all particles other than the lepton were clustered into two 
jets using the DURHAM algorithm. These two jets were required to contain at least 
three particles, at least one of which had to be charged. Additional preselection cuts arc 
listed in Table 4: these reject most events due to photon-photon collisions, some events 
for which there is no missing energy, and events whose topologies are far from those of 
WW events. 

After these preselection cuts, the final selection was made with an Iterative Discrim- 
inant Analysis (IDA) [27]. The standard IDA technique assumes that the signal and 
background distributions in the multi-observable space have different means but iden- 
tical shapes. IDA was extended to treat correctly cases when the distributions have 
different shapes. The input observables were transformed to make their distributions 
Gaussian. The IDA selection was made in three channels {qq^v^ qqev and qqrv). The 
training was performed on Monte Garlo samples: around 50k events each of four-fermion 

■^The isolation angle Biso is defined as the angle made to the closest charged particle with a momentum greater than 1 
GeV/c. 

■^defined as ^^^{^^• \pi\)/ \'Pi\> where 6i is the angle made by the momentum pi of the i"' particle in the jet with 
the total jet momentum 
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qqeu/qqixu 


qqru 


Transverse energy (GeV) 
Missing momentum (GeV/c) 
Visible energy /Etrain (%) 
Fitted W mass (GeV/c^) 


> 45 

> 10 

40<E/Et,,,„<110 

> 50 


> 40 

10<p^i,<80 
35<E/Et,,,„<100 

> 50 



Table 4: Semi-leptonic preselection cuts. Etrain indicates the centre-of-mass energy cho- 
sen for the training of the IDA (see text). The W mass is fitted through a constrained 
fit of the measured jets, lepton and unobserved neutrino imposing four-momentum con- 
servation and the equality of the two W masses in the event. 

charged and neutral current processes, and 100k qq{'y) events. The IDA was trained at 
three centre-of-mass energies: 189, 200 and 206 GeV. 

The observables used in the discriminants are summarised in Table 5, and can be 
grouped into four categories: 

• event observables such as the multiplicity, visible and transverse energies, aplanarity^ 
and the reconstructed effective centre-of-mass energy (v^) are useful to distinguish 
between semi-leptonic and other four-fermion events, and to remove qq{'y) events, 
particularly those in the radiative return peak; 

• observables concerned with the charged lepton (energy, isolation angle, tau jet mul- 
tiplicity) are useful in discriminating between events with primary leptons and those 
with leptons originating in quark decays or in other processes, such as photon con- 
version; 

• properties of the missing momentum (magnitude and polar angle) make use of the 
fact that the missing momentum from unseen initial state radiation (ISR in what 
follows) photons is directed along the beampipe; 

• the angles between the lepton and missing momentum, and the fitted W mass from 
a constrained kinematic fit (imposing four-momentum conservation and the equality 
of the two W masses in the event) are sensitive to the event topologies expected 
from semi-leptonic WW decays. 

The performance of the selection was measured using samples independent from those 
on which the training was performed. Events were selected with a cut on the output of 
the IDA, chosen to optimise the product of efficiency and purity for each channel. Events 
were first passed to the qq/iu selection; if they were not selected, they were passed to the 
qqeu; if they were still not selected, they were then finally passed to the qqru selection. 

Distributions of discriminants for the semi-leptonic event selection are shown in Fig- 
ure 3. 

5.2.2 Results for semi-leptonic final state 

The efficiency matrices and background contaminations for the semi-leptonic event 
selections were evaluated independently at the different centre-of-mass energies. The 
efficiencies differed by no more than 2% over the different energy points. The values 

*the aplanaxity is defined as I A3 , where A3 is the smallest eigenvalue of the sphericity tensor S"^ = , ' — ^ . The pi 

2^. iPiP 

are the 3-momenta of particles in the event, and a,j3 = 1,2,3 correspond to the x,y,z momentum components. 
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qqeu 


qqflu 


qqru 


total multiplicity 


X 


X 


X 


visible energy 


X 


X 


X 


transverse energy 


X 






aplanarity 






X 


\/?/^ 


X 


X 


X 


lepton energy 


X 


X 




Icpton isolation angle 


X 


X 


X 


tan jet charged multiplicity 






X 


tan jet momentum-weighted spread 






X 


magnitude of missing momentum 


X 


X 


X 


polar angle of missing momentum 


X 


X 


X 


lepton-missing momentum angle 




X 


X 


Fitted W mass 


X 


X 


X 



Table 5: Observables in semi-leptonic channels used in IDA analysis. 

at 200 GeV, with the number of events observed at each energy point are reported in 
Table 2. The final efficiencies for identifying qqeu, qqfiu and qqru events with any of the 
three semi-leptonic selections at 200 GeV were 79.6%, 89.6% and 63.9% respectively. The 
background in all the channels was composed of two- and four-fermion events in similar 
proportions. 

The number of events observed in the different lepton channels was found to be consis- 
tent with lepton universality. With this assumption, likelihood fits to the Poisson prob- 
ability of the expected number of events, where the contamination of the fully hadronic 
channel is considered as background, yield the cross-sections cr^w — '^ww x BR(WW —>■ 
qqliy) reported in Table 6. Systematic uncertainties were determined as detailed in Sec- 
tion 6. 

5.3 Fully leptonic final state 

5.3.1 Selection of fully leptonic final state event 

Events in which both W bosons decay into li/ are characterised by low multiplicity, a 
clean two-particle (or jet from r decay) topology with two energetic, acoUinear leptons 
of opposite charge, acoplanar with the beam and with large missing momentum and 
energy. The relevant backgrounds are di-leptons from e+e~ — > Z(7), Bhabha scattering, 
two-photon collisions and leptonic final states from Z-pair and single boson production. 

The selection was performed in three steps. First a leptonic preselection was made, 
followed by particle identification. Finally two Neural Networks were used to reject the 
remaining background. 

An initial set of cuts was applied to select a sample enriched in leptonic events. 
All particles in the event were clustered into "jets" using the LUCLUS algorithm 
{djoin — 6.5 GeV/c) and only events with two reconstructed jets, containing at least 
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(GeV) 


a^w = ^ww X BR(WW ^ 


qq\iy) (pb) 


161 


l.ntoil (stat) ±0.10 (syst) 


172 


5.68tlf3 (stat) ± 0.17 ( 


svst) 


183 


7.25 ± 0.46 (stat) ± 0.08 


fsvst) 


189 


6.89 ±0.26 (stat) ±0.08 


(syst) 


192 


6.88 ± 0.64 (stat) ± 0.08 


(syst) 


196 


7.50 ± 0.39 (stat) ± 0.08 


(syst) 


200 


7.82 ± 0.37 (stat) ± 0.08 


(syst) 


202 


7.71 ± 0.53 (stat) ± 0.08 


(syst) 


205 


7.45 ± 0.38 (stat) ± 0.08 


(syst) 


207 


6.95 ± 0.29 (stat) ± 0.08 


(syst) 



Table 6: Measured semi-leptonic cross-sections. 



one charged particle each, were retained. A total charged particle multiplicity between 2 
and 6 was required and at least one jet had to have only one charged particle. In order 
to reduce the background from two-photon collisions and radiative di-lepton events, the 
event acoplanarity, Qacop^ defined as the acollinearity of the two jet directions projected 
onto the plane perpendicular to the beam axis, had to be above 3°. In addition, the total 
momentum transverse to the beam direction, Pj, had to exceed 2% of the centre-of-mass 
energy ^/s. The associated energy deposited in the electromagnetic calorimeters for both 
leading particles (the ones with the largest momenta) was required to be less than 0.44- 
to reject Bhabha scattering. To reject radiative events further, the energy of the most 
energetic photon had to be less than 0.25 ■ ^/s and the angle in the plane perpendicular 
to the beam axis between the charged particles system and the most energetic photon 
was required to be less than 170°. Finally, the energy of the charged particles in each jet 
had to be greater than 0.04 • ^/s and the visible energy of the particles with | cos^ |< 0.9 
had to exceed 0.06 ■ a/s. 

In events passing this selection each particle was classed as /x, e or hadron. A lepton 
was identified as a cascade decay from W — > rv^ if the momentum was lower than 0.13--\/i. 

After the preselection and the channel identification, two Neural Networks based on 
the Multi Layer Perceptron (MLP) [28] package were built to reject the remaining back- 
ground. They consisted of one output layer from 13 input variables, normalised to lie in 
the region between zero and one. 

One of the two Neural Networks was tuned for the tVt-Xvx channels {X = e, /x, r) 
and the other for the remaining channels, given the different characteristics of the two 
samples. This was found to optimise the performance of the selection. 

The following variables were used in both Neural Networks: 

• the event acoplanarity; 

• the event acollinearity; 

• the larger of the associated energy from electromagnetic calorimetry of the two 
leptons; 

• the transverse momentum, pt] 

• the transverse energy; 
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• the angle in the plane perpendicular to the beam axis between the vector sum of the 
charged particle momenta and the most energetic photon; 

• the absolute value of the cosine of the polar angle of the missing momentum; 

• the energy in the calorimeters not associated to charged particles outside two cones 
of 20° around both leading charged particles; 

• the larger of the energies of the two jets coming from charged particles; 

• the total energy of neutral particles; 

• the larger of the invariant masses of the two jets; 

• the total visible energy; 

• the total energy of charged particles. 

Distributions of the Neural Network output variable for the fully leptonic event selection 
are shown in Figure 4. The cut applied on the Neural Network output was tuned in order 
to optimise the product of efficiency and purity. The cut was channel dependent, but did 
not depend on the centre-of-mass energy. 

5.3.2 Results for fully leptonic final state 

The efficiencies and backgrounds for ^/s = 200 GeV together with the numbers of 
selected events at each of the centre-of-mass energies are shown in Table 7. The overall 
Ivlv efficiency was 67.3% and the residual background from non-W and single-W events 
was 0.187 pb. The efficiencies were checked to be constant within 1-2% at the different 
energy points above 172 GeV. The number of observed events in each subchannel is 
consistent with the hypothesis of lepton universality. Table 8 presents the values, at each 
energy, of ~ '^ww x BR(WW \v\v) from maximum likelihood fits to the Poisson 
probability of the expected number of events. Systematic uncertainties were determined 
as detailed in Section 6. 

6 Systematic errors 

A large variety of systematic effects were taken into account in the cross-section deter- 
mination. They were due to imperfect modelling of the detector response or of underlying 
physics in the simulation, or to statistical uncertainties due to the finite size of the sim- 
ulation samples. 

Systematic errors were divided in three classes to facilitate their treatment in the com- 
bination of cross-section results from the four LEP experiments [29]: errors correlated 
between the experiments and between the different centre-of-mass energies (LCEC), typ- 
ically due to the use of the same models to describe physics effects, errors uncorrelated 
between the experiments but correlated between centre-of-mass energies (LUEC), which 
comprise detector-related effects, and errors correlated neither between experiments nor 
between energies (LUEU) which are mainly due to uncertainties on simulation sample 
statistics. Correlations between the different channels were also taken into account. 

Details of the determination of the systematic contributions are given in the following. 
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efficiencies for selected channels 


channel 


TiyTiy 












TfTV 


0.272 


0.087 


0.077 


0.004 


0.008 


0.006 


evTv 


0.068 


0.462 


0.005 


0.044 


0.049 


< 2 • 10-3 


\XVTV 


0.042 


0.004 


0.536 


< 2 ■ 10-3 


0.061 


0.053 


evev 


0.017 


0.151 


< 10-3 


0.471 


< 10-3 


< 2 ■ 10-3 




0.012 


0.041 


0.094 


< 2 ■ 10-3 


0.621 


< 10-3 




0.008 


< 10-3 


0.109 


< 2 ■ 10-3 


0.003 


0.677 


background (pb) 


0.030 


0.046 


0.023 


0.042 


0.019 


0.027 


(GeV) 


Selected events 


161 








2 






172 








8 






183 


4 


13 


10 


9 


9 


14 


189 


22 


56 


43 


25 


45 


38 


192 


4 


2 


14 


4 


8 


7 


196 


16 


29 


19 


12 


24 


7 


200 


12 


26 


28 


11 


27 


13 


202 


4 


20 


13 


4 


6 


9 


205 


14 


26 


22 


7 


24 


10 


207 


14 


40 


41 


16 


44 


16 


All 


891 



Table 7: Data for the cross-section measurement of the fully leptonic final state. The 
efficiency matrix and the background are the ones at 200 GeV. The upper hmits on the 
efficiencies are at 95% C.L. 

6.1 Estimation of systematic uncertainties 

6.1.1 Background cross-sections from four-fermion and two-fermion pro- 
cesses 

Theoretical uncertainties in the knowledge of four-fermion cross-sections largely de- 
pended on the process, varying from ±2% in ZZ processes to ±5% in single-boson con- 
tributions. Larger uncertainties arose when considering regions of the phase space domi- 
nated by 77 collisions [29]. 

The four-fermion background generation was performed independently for charged and 
neutral current processes. The latter ones were further divided with phase space cuts 
to isolate as much as possible the contribution of multiperipheral diagrams, where 77 
scattering dominates, and which were generated in a complementary sample [17]. For the 
selected background from charged current processes, mainly the single-W contribution, a 
relative uncertainty of ±5% was assigned. The same uncertainty was used for the non- 
77 neutral current background, with the exception of ZZ contributions where it became 
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(GeV) 


a^^% = dww X BR(WW ^ (pb) 


161 


0.30t[!:^^ (stat) ± 0.09 (syst) 


172 


l-03t[J:i (stat) ±0.09 (syst) 


183 


1.59 ± 0.26 (stat) ± 0.08 (syst) 


189 


1.86 ± 0.14 (stat) ± 0.04 (syst) 


192 


1.97 ± 0.37 (stat) ± 0.06 (syst) 


196 


1.87 ±0.21 (stat) ±0.05 (syst) 


200 


1.84 ± 0.20 (stat) ± 0.05 (syst) 


202 


1.81 ± 0.28 (stat) ± 0.05 (syst) 


205 


1.82 ± 0.20 (stat) ± 0.06 (syst) 


207 


1.82 ± 0.16 (stat) ± 0.06 (syst) 



Table 8: Measured fully leptonic cross-sections. 



±2%. For the remaining neutral current phase space a ±10% relative uncertainty was 
assigned. 

The theory errors used for the two-fermion background cross-sections were those re- 
ported in [30]. A further uncertainty of ±1% was conservatively added to account for 
ungenerated regions of the phase space in the matching between two-fermion and four- 
fermion processes, where radiative corrections to ee — > // could convert into an extra 
low-mass fermion pair. 

The variations of each background were taken as fully correlated in channel and energy 
and the resulting change in the measured cross-section was quoted as a systematic error. 

6.1.2 Modelling of four-jet background from qq 

In the fully hadronic channel the uncertainty on the background from two-fermion 
production leading to four-jet final states was estimated by comparing simulations with 
different hadronisation models, ARIADNE, PYTHIA and HERWIG. ARIADNE was 
chosen as default, because it provided the best description of the four-jet rates observed for 
the large data set at the Z peak [20] and also at LEP2 energies [31]. For the other models 
differences of (-6.0 ± 0.4)% (PYTHIA) and (±3.9 ± 0.6)% (HERWIG) were obtained for 
the background computation, and the largest difference of 6% was conservatively chosen 
as the systematic uncertainty. This estimate was confirmed by fitting all data to signal 
and background with a free parameter to scale the two-fermion background (ARIADNE), 
for which a change of (4 ± 3)% was obtained. 

6.1.3 Fragmentation modelling 

Modelling of the fragmentation in hadronic events could have an impact both on the 
selection efficiency and on the background level estimated from the Monte Carlo. These 
effects were evaluated by comparing the performance of the selection algorithms on signal 

and background samples generated with different hadronisation models. 

The PYTHIA hadronisation model, which best described the two-jet fragmentation 
at the Z peak, was taken as a reference to evaluate efficiencies for all channels and 
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the background for the semi-leptonic channel. It was compared to ARIADNE [12] and 
HERWIG [13] and the largest deviations were considered to estimate the systematic errors 
on the cross-sections. 

The largest variation in the signal efficiency was 0.58% in the fully hadronic channel 
and 1% in the semi-leptonic one, while the difference in the background level was found 
to be 3.2% in the semi-leptonic channel. For the hadronic channel this uncertainty is 
included in the four-jet modelling. 

6.1.4 Final state interactions 

At LEP2 energies the decay distance between the W bosons was smaller than the 
hadronisation scale or the typical radius where Bose-Einstein effects occur. Therefore 
gluon exchange between quarks from different Ws (known as Colour Reconnection) or 
Bose-Einstein Correlations between pions were to be expected. These so-called Final 
State Interactions (FSI) between the decay products of the two different Ws could affect 
the reconstruction of fully hadronic events only, and their modelling could have an impact 
on the determination of the selection efficiency. 

The effects of Colour Reconnections and Bose-Einstein Correlations were estimated by 
evaluating the selection efficiency on simulation samples where FSI were modelled using 
the SKI algorithm [32] with reconnection probability of 30% for Colour Reconnections and 
the LUBOEI algorithm [33] for Bose-Einstein Correlations. The full difference between 
the presence and the absence of these effects was taken as an indication of the systematic 
error, corresponding to a variation in efficiency of -0.3% and 0.2%, respectively. 

6.1.5 Radiative corrections 

The correct simulation of radiation in WW production and decay could be relevant 
for efficiency determination. Since the LEP2 Monte Carlo workshop [1], generators 
with a more precise calculation of 0{a) electroweak corrections to CC03 became avail- 
able [34,16]. The theoretical uncertainty on the total cross-section was reduced by almost 
a factor 4 down to a level of 0.5%, with a change of central value by almost 2% with re- 
spect to Gentle2.0 [35], run with parameter settings as described in [36]. Moreover, it was 
shown that the more correct description of electroweak corrections also had very impor- 
tant effects on differential distributions [37], making its inclusion essential for efficiency 
determinations . 

The YFSWW [16] program was used in the generation for reweighting events accord- 
ing to the procedure described in [17]. In order to estimate the effect that theoretical 

uncertainty in the description of the radiation had on the selection efficiency, the re- 
sults obtained in the DELPHI setup were compared with a simulation making use of 
RacoonWW [34] on the testbench process udiiu^. The two programs differed in many 
respects in the treatment of QED radiation: YFSWW was a e+e~ W+W" — > 4f 
generator (CC03 diagrams only) with 0{a) factorisable electroweak corrections and non- 
factorisable corrections implemented via the so-called Khoze-Chapovsky ansatz (KC) [38]. 
It included ISR in leading logarithm approximation (LL) 0{a^) via YFS exponentia- 
tion [39] and final state radiation (FSR) LL 0{a^) via PHOTOS [40]. On the other 
hand, RacoonWW was a 4/ generator implementing 0{a) in Double Pole Approxima- 
tion (DPA) rigorously, extended to 0{a^) for collinear ISR via structure functions. It 
also included real corrections with the exact e+e~ — > 4/7 matrix elements of the CCll [2] 
class. 
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The systematic uncertainty on the cross-section was estimated by comparing results 
from the two programs: no significant difference was found in the efficiencies, therefore the 
statistical error on this difference was taken as a conservative estimate of the systematic 
contribution. 

6.1.6 Luminosity determination 

The luminosity was determined from a measurement of Bhabha scattering, which was 
theoretically known to high accuracy. The measurement made use of coincidences in the 
very forward electromagnetic calorimeters and was affected by the experimental error 
on the acceptance (±0.5%). The residual theoretical uncertainty on the cross-section 
estimate was ±0.12% [41]. These uncertainties were propagated to a systematic error on 
cww- 

The statistical error on the Bhabha cross-section was included in the statistical error 
of the W-pair cross-section. 

6.1.7 Detector effects 

A non-perfect reproduction of track reconstruction and lepton identification efficien- 
cies in the simulation could induce a systematic error on the signal efficiency and the 
background rejection. This was particularly relevant for semi-leptonic and fully leptonic 
channels. 

These effects were evaluated from a comparison of the simulation with data on high 
statistics samples of clean two-lepton and two-jet events, which were collected at a centre- 
of-mass energy of 91.2 GeV. The data were taken with the same detector and trigger 
configuration and analysed with the same reconstruction software as the high energy 
data. High energy data samples were also used for the comparison, selecting events with 
a clean di-lepton and 2- or 3-jet topologies. 

From data-simulation comparisons, corrections were deduced for the energies and polar 
angles of jets, and for the number of charged particle tracks in the endcap region. These 
led to corrections of the signal efficiency in the hadronic channel of ±0.14% (jets) and 
±0.30% (endcap tracks), for which systematic uncertainties of ±0.14% and ±0.10% were 
assigned. 

For lepton identification, the angular averaged difference in the identification efficiency 
between data and simulation of 0.3% for muons and 1% for electrons was found. The 
systematic error on the WW cross-section and W branching fractions was determined by 
randomly changing the lepton identification in the signal and background Monte Carlo 
according to the discrepancies found. The changes were correlated between the different 
channels. The effect on the total cross-section was small because most of the events lost 
as electrons or muons were recovered by the tau selections. 

Other possible discrepancies between data and Monte Carlo in the variables used 
as inputs for the Neural Network in the hadronic channel were taken into account by 
smearing these variables in the simulation by their experimental resolution. 

In total a systematic uncertainty of ±0.55% on the signal efficiency and ±1.7% on the 
total background was estimated for the hadronic channel. For the semi-leptonic channels 
the uncertainties on the signal efficiencies were ±0.5% for the qqev and qq^v channels, and 
±1% for the qqrv. The variations of the background ranged from 0.6% to 6% according 
to the channel. The systematic uncertainties in the fully leptonic channel were ±1.5% 
for the efficiency and 7% for the background. 
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6.1.8 Detector inefficiencies in a specific period 

As mentioned in Section 2, one of the TPC sectors, covering 1/12 of the acceptance, 
was not operational during the last period of the high energy data taking in the year 
2000. These data were analysed separately and then combined with the results from the 
previous period. The performance of the analyses and the cross-section values were found 
to be compatible within statistical errors. Additional systematic effects were estimated 
by comparing data collected at the Z peak during the period with the TPC sector 6 
off with simulation samples produced with the same detector conditions. Both hadronic 
and leptonic Z decays were used. The impact on the WW cross-section analysis was 
conservatively evaluated as an uncertainty on the selection efficiency of ±0.5% in the fully 
hadronic channel and of 1% in the other channels, which was added to the systematic 
error. No effects were found on the background level. 

6.1.9 Monte Carlo statistics 

The uncertainties due to limited statistics of the Monte Carlo samples were at most 
±0.2% on the signal efficiency and ±2.0% on the background level for the hadronic 
channels at energies above 172 GeV. The error on backgrounds in the fully leptonic 
channels were up to 10 times larger because of the small numbers of accepted events. 



Source 


(pb) 


^ww (pb) 


^l^i'^ (pb) 


LC 


EC 


Four-jet modelling 


±0.051 


±0.014 




Y 


Y 


Background cross-sections 


±0.009 


±0.016 


±0.006 


Y 


Y 


Fragmentation 


±0.045 


±0.038 




Y 


Y 


Final state interactions 


±0.025 






Y 


Y 


Radiative corrections 


±0.008 


±0.008 


±0.002 


Y 


Y 


Luminosity (theor) 


±0.011 


±0.010 


±0.002 


Y 


Y 


Luminosity (cxp) 


±0.045 


±0.043 


±0.011 




Y 


Detector effects 


±0.045 


±0.053 


±0.033 




Y 


Monte Carlo statistics 


±0.005 


±0.014 


±0.033 







Table 9: Breakdown of systematic errors on the partial WW cross-sections at v^=200 
GeV and their classification according to correlations between LEP experiments (LC) 
and between centre-of-mass energies (EC). 

Table 9 presents the breakdown of the systematics per channel determined in the way 
described above for 200 GeV. The classification of the systematics according to their 
correlations in energy or experiments is also shown. 

7 Differential cross-section 

Given the high statistics available at LEP2, it is interesting to provide a measurement 
of differential cross-sections. Particularly relevant are the W polar angle distributions, 
from which triple gauge coupling limits can be extracted. In what follows the determi- 
nation of d(Jww/dcos6'w- is discussed. 
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The presented differential cross-section refer to CC03 qqev and qq/j.i' final states, 
since they give a clean W charge assignment, and has therefore to be understood 
as d[(Tww(BReiy+BR^,^)]/dcos^w-, where BRi^, is the branching ratio of the decay 
WW ^ qqli^. 

The IDA analysis described in section 5.2.1 was used for the event selection; the W 
flight direction was determined via a constrained flt of the event which imposed four- 
momentum conservation and equality of the hadronic and leptonic invariant masses. In 
order to distinguish between radiation from Ws or from the final state fermions, the 
photon-to-fermion recombination scheme followed the CAL05 definition adopted in the 
LEP2 Monte Carlo workshop [1]. To match the detector acceptance best, an additional 
restriction requiring the charged lepton to be more than 20° away from the beam direction 
was introduced. For the qqeu case this additional cut strongly suppressed the contribution 
of single-W diagrams in the signal definition. The impact of possible systematic effects 
due to charge misidentification for leptons within the accepted angular region was found 
to be negligible from studies on dilepton events. To optimise the statistics in each bin, 
four bins in energy and ten bins in polar angle were chosen. Table 10 reports the energy 
binning, together with the corresponding luminosity and luminosity weighted centre- 
of-mass energy. The angular binning was chosen to be significantly bigger than the 
resolution on cos6'w-, estimated to be about 0.06 from Monte Carlo, in order to minimise 
bin migration of data. The migration matrix, expressing the probability that an event 
selected in a certain bin was generated in another one, is reported in Appendix A. No 
corrections for bin migration were applied to the presented results. All the conventions 
used for presenting the DELPHI result, including signal and bin definitions, are the ones 
agreed for the LEP combination. 



Bin (GeV) 


/^-weighted ^/s (GeV) 


C (pb-i) 


180-184 


182.65 


51.63 


184-194 


189.03 


178.32 


194-204 


198.46 


193.52 


204-210 


205.91 


198.59 



Table 10: Energies and luminosities in the bins for the differential cross-section measure- 
ment. 

The results at the energies in Table 10 are reported in Figure 5, where data points are 

superimposed on the expected distributions from WPHACT and YFSWW. The data are 
in agreement with the expectations in all energy ranges. A systematic deficit of data in 
the highest cos6'w- bin at energies above 184 GeV is observed. As a crosscheck, it was 
verified that the shape of the angular distributions, for electrons and muons separately, 
and for positively and negatively charged leptons, were compatible within errors. 

The detailed list of results, in terms of measured cross-sections, statistical and sys- 
tematic errors per bin is reported in Appendix A. 
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8 Determination of the W branching fractions and 
of the total WW production cross-section 

The total cross-section for WW production and the W leptonic branching fractions 
were obtained from a likelihood fit based on the probabilities of finding the observed num- 
ber of events in each final state. The input numbers in the form given in Tables 2 and 7 
were used, except for the fully hadronic final state where, for energies above 172 GeV, 
the binned distribution of the neural network output was used. 

A fit without the assumption of lepton universality, requiring the unitarity of the 
branching ratios (i.e. BRgi^ + BR^,y -|- BKru + BR^-/ = 1), was performed. The results 
for all the data above 172 GeV arc shown in Table 11. The analyses on the low statistics 
samples at 161 and 172 GeV used inclusive lepton identification. The correlation matrix 
is also reported, where both statistical and systematic contributions are included. Since 
the lepton branching ratios are in agreement, a second fit assuming lepton universality 
was performed in order to extract BR^-/. The result for this second fit are also reported 
in Table 11. In Figure 6 the results are compared to the Standard Model predictions. All 
results are consistent with the expectations and with lepton universality. 



channel 


branching fraction 


stat. error 


syst. error (LU) syst. error (LC) 


W ^eu 

W ^ I2U 
W ^TV 


0.1055 
0.1065 
0.1146 


0.0031 
0.0026 
0.0039 


0.0013 0.0005 
0.0006 0.0005 
0.0017 0.0009 




Correlations 


W ^ ez/ W ^ fiu W ^ TV 




W ^ei/ 
W ^ l^u 

W^Tl/ 


1.00 0.03 -0.34 
0.03 1.00 -0.17 
-0.34 -0.17 1.00 





assuming lepton universality 


channel 


branching fraction 


stat. error 


syst. error (LU) syst. error (LC) 


W — > hadrons 


0.6745 


0.0041 


0.0020 0.0014 



Table 11: W branching fractions from data above 172 GeV and correlation matrix for 
the leptonic branching fractions. 



Within the Standard Model, the branching fractions of the W boson depend on the six 
matrix elements |Vqq'| of the Cabibbo-Kobayashi-Maskawa (CKM) quark mixing matrix 
not involving the top quark. In terms of these matrix elements, the leptonic branching 
fraction of the W boson BR;,^ is given by 



1 



BRiu 



3U + 



1 + 



".(M^v) 



TT 
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where ^^(M^) is the strong couphng constant. Taking ^^(M^) = 0.119 ±0.002 [42], the 
measured leptonic branching fraction of the W yields 

J2 |Vij|^ = 1.996 ± 0.043 (BRi^) ± 0.002 (as), 

i = (u, c), 
3 = (d,a,V) 

where the first error is due to the uncertainty on the branching fraction measurement 
and the second to the uncertainty on a^- Using the experimental knowledge [42] of the 
sum |VudP + |VusP + |VubP + |VcdP + |VcbP = 1.0476 ± 0.0074, the results of Table 11 
can be interpreted as a measurement of |Vcs| which is the least well determined of these 
matrix elements: 

|Vcs| = 0.973 ± 0.019(stat) ± 0.012(syst) 

where the uncertainties on the SM parameters are included in the systematic error. 

If the SM values of the W branching ratio [42] are assumed, the fitting procedure can 
be repeated with the total WW production cross-sections as the only free parameter. 
The results obtained from all the data samples are reported in Table 12. In this table 
the breakdown of the systematic contributions into the correlation categories defined in 
Table 9 is also shown. Correlations between the different channels were also taken into 
account. 



x/i (GeV) 


aww (pb) 








'syst.) 


(pb) 








LCEC 


LUEU 


LUEC 


161 


3.6llg|^ (stat) ± 0.19 (syst) 





039 


0.182 


0.037 


172 


11.37;J:^^ (stat) ±0.32 ( 


syst) 





115 


0.288 


0.095 


183 


16.07 ±0.68 (stat) ±0.16 


(syst) 





100 


0.044 


0.126 


189 


16.09 ±0.39 (stat) ±0.16 


(syst) 





094 


0.028 


0.127 


192 


16.64 ±0.99 (stat) ±0.17 


(syst) 





103 


0.040 


0.131 


196 


17.04 ±0.58 (stat) ±0.17 


(syst) 





100 


0.033 


0.133 


200 


17.39 ±0.55 (stat) ±0.17 


(syst) 





102 


0.029 


0.135 


202 


17.37 ±0.79 (stat) ±0.17 


(syst) 





107 


0.030 


0.136 


205 


17.56 ±0.57 (stat) ±0.17 


(syst) 





107 


0.039 


0.136 


207 


16.35 ± 0.44 (stat) ± 0.17 


(syst) 





103 


0.051 


0.128 



Table 12: Measured total WW cross-sections. A breakdown of the systematic uncertain- 
ties in the different contributions, correlated between LEP experiments and/or between 
different centre-of-mass energies, is shown, according to the classification described in the 
text (see Section 6 for details). 

A comparison of the results with the most recent calculations in the Double Pole 
Approximation from RacoonWW and YFSWW is shown in Figure 7. As DPA compu- 
tations are not rehable close to the WW threshold, the predictions below 168 GeV were 
obtained running those programs in the Improved Born Approximation (IBA), which 
only accounts for initial state radiation and Coulomb corrections. The shaded region 
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represents the theoretical uncertainty of the calculations and was obtained by an ana- 
lytic parametrisation of the relative uncertainty given in [1]. This led to an accuracy 
on the theoretical curves of about ±0.7% at 168 GeV and of ±0.4% at 200 GeV. The 
uncertainties from RacoonWW and from YFSWW were merged into a single error band. 
The measurements are in very good agreement with the Standard Model expectations. 

9 Determination of T^ww 

The cross-section measurements for different energies can be combined into a single 
value to quantify the overall agreement with theoretical predictions. T^-ww is defined 

as the ratio between the experimentally determined cross-section and the theoretical 
expectations. This procedure was used to compare the measurements at the eight ener- 
gies between 183 and 207 GeV to the predictions of Gentle, KoralW, YFSWW and 
RacoonWW. The measurements at 161 and 172 GeV were not used because of the high 
sensitivity of the cross-section to the value of the W mass at these energies. 

For each calculation, the cross-sections were converted into ratios by dividing them by 
the corresponding theoretical predictions, and combined taking into account the energy 
correlation of the systematics. 

The theoretical errors on the calculations, due to both physical and technical precision 
of the generators used, were also propagated to the ratios. 



Theoretical prediction 


7?.ww 


Gentle2.0 
KoralW 
YFSWW 
RacoonWW 


0.974 ± 0.015 (exp) ± 0.019(theo) 
0.979 ± 0.015 (exp) ± O.OlO(theo) 
0.999 ± 0.015 (exp) ± 0.005(theo) 
1.001 ± 0.015 (exp) ± 0.005(theo) 



Table 13: T^ww values from the combination of all the data, using different theoretical 
calculations. 

The values of 7?-ww at the various centre-of-mass energies are presented in Figure 8 
for the RacoonWW calculation. The band on the figure represents the theoretical error 

on the prediction, where its dependence on energy is ignored for simplicity. The T^ww 
values from the combination of all the data using different theoretical calculations are 
shown in Table 13. Both statistical and systematic contributions are indicated. It is worth 
noting how the estimated theory error (taken from [1], page 34) decreases from Gentle to 
RacoonWW and that the final experimental precision on T^ww with the DELPHI data at 
LEP2 approaches ±1.5%, to be compared with ±0.5% from the best theory computations. 
The data favour the more complete inclusion of radiative corrections in the calculation 
at the level of 1.5 cr. 

10 Conclusions 

The WW production cross-section from e+e" annihilations has been measured at ten 
centre-of-mass energies between 161 and 209 GeV with the DELPHI experiment at LEP. 
The data correspond to a total integrated luminosity of about 670 pb~^. 
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The results are in agreement with the predictions of the most recent CC03 cross-section 
calculations and test them to an accuracy of about 1.5%. 

Differential distributions in the polar angle of the reconstructed W are also measured 
in the semileptonic channels. 

The W branching fractions are measured with an uncertainty of less than ±4% for 
individual leptons and about ±0.7% for hadrons. They also agree with the Standard 
Model expectation. Prom the leptonic branching fraction a precise determination of \Vcs\ 
is derived. 
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Appendix A: Results on W polar angle measurements 

The details of the results for the measurement of the W~ polar angle differential 
cross-section are reported in Table 14. For convenience the integrated luminosities and 
the lumosity-weighted centre-of-mass energies in each of the energy bins defined in the 
text are also provided. 

Table 15 presents the signal event migration matrix determined from the Monte Carlo. 
The Standard Model W polar angle differential distribution is therefore imphcitly as- 
sumed. 
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180-184 GeV 


184-194 GeV 


194-204 GeV 


204-210 GeV 


C (pb-1) 


51.63 


178.32 


193.52 


198.59 


weighted ^/s (GeV) 


182.65 


189.03 


198.46 


205.91 


ai [-1,-0.8) 


0.715 


0.865 


0.600 


0.275 


(5(Tf "*'(meas) 


(Jaf (exp) 


0.320 


0.320 


0.180 


0.165 


0.155 


0.150 


0.120 


0.145 


5(ji^**'(back) 


5(rr*'(eff) 


0.045 


0.020 


0.040 


0.020 


0.025 


0.015 


0.020 


0.015 


(72 [-0.8,-0.6) 


0.795 


0.760 


0.675 


0.590 


5(T|*°'*-(meas) 




0.315 


0.315 


0.170 


0.170 


0.160 


0.160 


0.145 


0.150 


5(T2^**'(back) 


5(7r*-(eff) 


0.025 


0.025 


0.020 


0.020 


0.015 


0.020 


0.015 


0.020 


CT3 [-0.6,-0.4) 


1.175 


0.990 


1.510 


0.575 


5(T3*"*-(meas) 




0.380 


0.350 


0.185 


0.180 


0.215 


0.170 


0.140 


0.160 


5cr3^'**'(back) 


5ar*-(eff) 


0.020 


0.035 


0.020 


0.035 


0.015 


0.030 


0.010 


0.025 


(74 [-0.4,-0.2) 


1.365 


0.930 


1.150 


0.930 


5o'4*"*-(meas) 


^^stat. (^Qxp) 


0.400 


0.370 


0.180 


0.200 


0.190 


0.180 


0.170 


0.175 


(5cr4^**'(back) 




0.015 


0.035 


0.015 


0.035 


0.015 


0.035 


0.010 


0.035 


[-0.2,0) 


1.350 


1.330 


1.055 


1.000 


(5(Tg*"*- (meas) 


§iyStat. (cxp) 


0.400 


0.405 


0.215 


0.215 


0.185 


0.200 


0.175 


0.195 


(5c7f'*-(back) 


5(7f^*-(eff) 


0.015 


0.040 


0.015 


0.040 


0.015 


0.035 


0.015 


0.035 


f76 [0,0.2) 


1.745 


1.460 


1.635 


1.190 




5(7g*"*-(exp) 


0.450 


0.450 


0.225 


0.240 


0.225 


0.230 


0.195 


0.220 


(^(7r*(back) 




0.025 


0.085 


0.020 


0.085 


0.015 


0.085 


0.010 


0.085 


(77 [0.2,0.4) 


1.995 


1.675 


2.115 


2.120 


(5 erf"*- (meas) 


(5(7f"*-(exp) 


0.485 


0.505 


0.240 


0.270 


0.255 


0.260 


0.255 


0.250 


(5(7r*(back) 


(5(7f^*-(eff) 


0.015 


0.050 


0.015 


0.050 


0.010 


0.045 


0.010 


0.045 


(78 [0.4,0.6) 


2.150 


2.630 


3.175 


2.655 


5(7;;"''-(iiieas) 


(ya^'""-(exp) 


0.510 


0.580 


0.300 


0.320 


0.320 


0.310 


0.290 


0.300 


(5(7r*(back) 


(5(7r*-(efF) 


0.015 


0.065 


0.015 


0.060 


0.015 


0.055 


0.010 


0.055 


(79 [0.6,0.8) 


4.750 


4.635 


4.470 


4.585 


(5(79*"*- (meas) 


(5(7|*"*-(exp) 


0.775 


0.695 


0.405 


0.385 


0.385 


0.380 


0.385 


0.380 


(5(7r*-(back) 


(5(7^^'*(eff) 


0.030 


0.095 


0.025 


0.100 


0.025 


0.105 


0.020 


0.110 


(710 [0.8,1] 


6.040 


5.400 


7.140 


7.290 




5c7iQ"*-(exp) 


0.895 


0.850 


0.455 


0.490 


0.500 


0.505 


0.505 


0.520 


5(7^0**' (back) 


<5(7r-(eff) 


0.035 


0.075 


0.035 


0.085 


0.030 


0.100 


0.030 


0.110 



Table 14: Differential cross-sections in the 10 angular bins for the four energy intervals 
(see descriptions in Section 7). (7j indicates the average of d[(7-ww(BRei/-l-BR^,y)]/dcos^w- 
in the i-th bin of cos6'w- with width 0.2. The limits of the bins are also reported in the 
table. The values, in each bin, of the measured and expected statistical error and of 
the systematic errors due to the background and to the efficiencies are reported as well. 
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gen 1 


gen 2 


gen 3 


gen 4 


gen 5 


gen 6 


gen 7 


gen 8 


gen 9 


gen 10 


sel 1 


71.8 


13.4 


2.1 


1.3 


0.5 


0.4 


0.1 


0.5 


1.8 


8.0 


sel 2 


11.4 


61.7 


15.4 


2.7 


2.1 


0.8 


0.6 


0.7 


2.4 


2.2 


sel 3 


1.5 


12.2 


57.2 


16.7 


4.2 


1.5 


1.4 


1.8 


2.1 


1.3 


sel 4 


0.8 


1.9 


11.7 


58.9 


16.8 


4.1 


2.6 


1.5 


0.9 


0.8 


sel 5 


0.2 


0.5 


1.7 


12.8 


56.4 


18.9 


5.0 


2.4 


1.1 


1.0 


sel 6 


0.2 


0.2 


0.4 


1.8 


12.9 


58.8 


18.7 


3.9 


1.9 


1.3 


sel 7 


0.1 


0.1 


0.4 


0.6 


1.7 


12.0 


59.8 


19.9 


4.0 


1.5 


sel 8 


0.0 


0.1 


0.1 


0.3 


0.5 


1.7 


11.6 


63.3 


19.4 


3.2 


sel 9 


0.1 


0.0 


0.1 


0.1 


0.2 


0.4 


0.9 


9.9 


70.7 


17.6 


sel 10 


0.1 


0.0 


0.0 


0.0 


0.1 


0.1 


0.2 


0.7 


7.4 


91.4 



Table 15: Polar angle migration matrix at 200 GeV for selected events of the sig- 
nal, according to the definition described in the text. The matrix is defined as 
Mij = Nseii;genj/NseH, where Nseii;genj IS the number of events selected in bin i and 
generated in bin j and Nggn is the total number of events selected in bin i. The numbers 
are all expressed in percent. By construction the rows sum up to 100%. The relative 
errors on the numbers on the diagonal are below 2%, whereas outside the diagonal they 
reach 7%. 
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Figure 1: CC03 diagrams. 




neural net output neural net output 

Figure 2: Distribution of the Neural Network output variable for four-jet events at the 
centre-of-mass energies of 189 GeV and 207 GeV. The period of bad TPC operating 
conditions is excluded from the plot. The points show the data and the histograms are 
the predicted distributions for signal and background. The arrows indicate the cut value 
applied for the selection of events. 



28 




-15 -12.5 -10 -7.5 -S -2.5 2.5 5 7.5 10 -10 -7.5 -5 -2.5 2.5 5 7.5 10 -IS -10 -5 5 10 15 

2nd qqev discriminant 2nd qq^v discriminant 2nd qqxv discriminant 



Figure 3: Distribution of discriminants for the semi-leptonic selection. The plots refer 
to the qqeu selection at 189 GeV, the ggfiu selection at 200 GeV and the qgru selection 
at 207 GeV. The arrows indicate the cut value applied for the selection of events. The 
points show the data and the histograms are the predicted distributions for signal and 
background. 
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Figure 4: Distribution of the output variable of the two types of Neural Network used 
for the selection of fully-leptonic events. The points show the data and the histograms 
are the predicted distributions for signal and background. 
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Figure 5: W~ polar angle differential cross-sections from the qqeu and gg/iz/ channels 
(detailed definition in the text). The measurements in the four energy bins defined 
in the text (points) are compared with the expectations (histograms). The systematic 
contributions are included in the error bars. 
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DELPHI W decay Branching Ratios 
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Figure 6: W decay branching ratios (without and with the assumption of lepton uni- 
versality) measured with the DELPHI data, in comparison with the Standard Model 
expectations [42] . Both the statistical and the total errors are indicated in the error bars. 
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Figure 7: Measurements of the WW cross-section compared with the Standard Model 
prediction given by the YFSWW [16] and RacoonWW [34] programs. The shaded band 
represents the uncertainty on the theoretical calculations. 



32 



DELPHI 



183 GeV 
189 GeV 
192 GeV 
196 GeV 
200 GeV 
202 GeV 
205 GeV 
207 GeV 

Combined 



I 



1.045 + 0.046 
0.990 ±0.026 
1.007 + 0.061 
1.014 ±0.037 
1.024 ±0.034 
1.020 ±0.048 
1.028 ±0.035 
0.957 ±0.028 



1.001 ±0.015 

^yndf = 6.54/7 



0.9 



1. 



1.1 



Figure 8: Ratios between measured and predicted CC03 cross-sections with the DEL- 
PHI data. The error bars indicate the combined statistical and systematic uncertain- 
ties. The energy-combined values are indicated at the bottom. The plot refers to the 
RacoonWW predictions, where a 0.5% theory error is indicated as a band. 



